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ABSTRACT. The separability of the massive Dirac equation in the non-extreme Kerr geometry in 
horizon-penetrating advanced Eddington-Finkelstein-type coordinates is shown. To this end, Kerr 
geometry is described by a Carter tetrad and the Dirac spinors and matrices are given in a chiral 
Newman-Penrose dyad representation. Applying Chandrasekhar’s mode ansatz, the Dirac equation 
is separated into systems of radial and angular ordinary differential equations. Asymptotic radial 
solutions at infinity, the event horizon, and the Cauchy horizon are explicitly derived. Their decay is 
analyzed by means of error estimates. Moreover, the eigenfunctions and eigenvalues of the angular 
system are discussed. Finally, as an application, the scattering of Dirac waves by the gravitational 
field of a Kerr black hole is studied. This work provides the basis for a Hamiltonian formulation 
of the massive Dirac equation in Kerr geometry in horizon-penetrating coordinates and for the 
construction of a functional analytic integral representation of the Dirac propagator. 
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I. INTRODUCTION 

Over the last five decades, the dynamics of relativistic spin-i fermions (Dirac waves) in Kerr geometry was studied 
extensively by employing different approaches. The probably most established approach is Chandrasekhar’s mode 
analysis [6-8], where the massive Dirac equation is separated by means of time and azimuthal angle modes and 
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rewritten in terms of radial wave equations and coupled angular ordinary differential equations (ODEs). Within 
this framework, many physical processes like the emission and absorption of Dirac waves by rotating black holes 
or black hole stability under fermionic field perturbations were investigated [3, 27, 29-31]. The dynamics of 
Dirac waves was also analyzed in the framework of scattering theory [1, 5, 20, 22], More recently, a functional 
analytic integral representation of the Dirac propagator was derived using the Hamiltonian formulation of the 
massive Dirac equation [15-17]. 

The basis of the mode analysis approach is Chandrasekhar’s famous discovery that the massive Dirac equation 
in Kerr geometry expressed in Boyer-Lindquist coordinates is separable, which was worked out in his original 
article from 1976 [6] and led to a major breakthrough in the field. At that time, this remarkable result came a 
bit as a surprise because the Dirac system of coupled Hrst-order partial differential equations (PDEs) was not 
expected to be separable in Kerr geometry. Despite the tremendous impact of this discovery, the validity of the 
solutions is naturally restricted to those regions of Kerr geometry where the Boyer-Lindquist coordinates are 
well-dehned. As they have singularities at the event and the Cauchy horizon, respectively, the dynamics of Dirac 
waves near and across these horizons cannot be described properly. In this article, we resolve this problem by 
applying a specific analytic extension of the Boyer-Lindquist coordinates that is regular at the horizons, thus, 
it covers both the exterior and interior black hole regions and allows for well-defined transitions of Dirac waves 
across the horizons. These so-called advanced Eddington-Einkelstein-type coordinates also feature a coordinate 
time required for the Hamiltonian formulation of the Dirac equation and the corresponding Cauchy problem. We 
point out that since the transformation to these coordinates is singular at the horizons and, hence, non-trivial, 
a careful analysis is necessary. Furthermore, as the mixing of the time and the azimuthal angle variable arising 
in this transformation leads to a symmetry breaking of structures inherent to Boyer-Lindquist coordinates, it 
is a priori not clear that the separation of variables property is conserved in advanced Eddington-Einkelstein- 
type coordinates. In the recent article [II], a similar coordinate system is employed. This work concentrates 
on the physical aspects of the massive Dirac equation in horizon-penetrating coordinates. More precisely, using 
numerical methods, the energy spectrum and decay rates of bound states are computed. Here, on the other hand, 
the focus is on the mathematical aspects, i.e., in the framework of the Newman-Penrose formalism, a detailed 
analysis of the asymptotics of Dirac waves at infinity, the event horizon, and the Cauchy horizon is given. 
This includes error estimates and, thus, the proper study of their decay. (For a different horizon-penetrating 
coordinate system see, e.g., [12].) 

In more detail, we perform the mode analysis of massive Dirac waves in the non-extreme Kerr geometry in 
horizon-penetrating coordinates as follows. First, in Section H, we describe Kerr geometry in the Newman- 
Penrose formalism by a regular Carter tetrad in advanced Eddington-Einkelstein-type coordinates and calculate 
the associated spin coefficients. Secondly, in Section HI A, we formulate the massive Dirac equation in the chiral 
representation with a Newman-Penrose dyad basis for the spinor space. General overviews of the Newman- 
Penrose formalism and of the general relativistic Dirac equation are given in Appendices A and B. Note that the 
Newman-Penrose formalism is well suited for the analysis of radiative transport in curved spacetimes, especially 
Dirac wave propagation in Kerr geometry because it can be chosen to reflect symmetries - or be adapted 
to certain aspects - of the underlying spacetime, which subsequently leads to a reduction in the number of 
conditional equations and to simplified expressions for geometric quantities. Next, considering a factorization of 
the Dirac waves in time and azimuthal angle modes, we show the separability of the massive Dirac equation in 
advanced Eddington-Einkelstein-type coordinates into systems of radial and angular ODEs. Asymptotic radial 
solutions at infinity, the event horizon, and the Cauchy horizon are determined in Sections HIB-HID. Moreover, 
error estimates that demonstrate suitable decay of these solutions are given. The angular system yields the 
usual Chandrasekhar-Page equation. The corresponding set of smooth eigenfunctions and the discrete, non¬ 
degenerate spectrum of eigenvalues are briefly, however appropriately, discussed in Section HIE. Finally, in 
Section IV, we apply the radial asymptotics at inhnity and at the event horizon to the physical problem of 
scattering of Dirac waves by the gravitational field of a Kerr black hole. To this end, we evaluate the net 
current of Dirac waves at space-like inhnity and at the event horizon. The resulting conservation law for the 
rehexion and transmission coefficients obtained for horizon-penetrating coordinates is in agreement with that 
for Boyer-Lindquist coordinates. Therefore, on the one hand, the conserved net current stays positive across 
the event horizon and, on the other hand, superradiance does not occur. This article provides the basis for a 
Hamiltonian formulation of the massive Dirac equation in the non-extreme Kerr geometry in horizon-penetrating 
coordinates. Within this framework, one can construct a functional analytic integral representation of the Dirac 
propagator that can be applied to Dirac waves with compact support in both the exterior and interior black hole 
regions. These problems are worked out in [18]. We point out that by using horizon-penetrating coordinates in 
the construction of the propagator, gluing techniques to connect the exterior and interior regions of the black 
hole, as employed in the Boyer-Lindquist case, are no longer required. 



II. NEWMAN-PENROSE REPRESENTATION OF KERR GEOMETRY IN 
HORIZON-PENETRATING ADVANCED EDDINGTON-FINKELSTEIN-TYPE COORDINATES 


The non-extreme Kerr geometry [23] is a connected, orientable and time-orientable, smooth, asymptotically flat 
Lorentzian 4-manifold (DJt, g) with topology 8“^ x that consists of a differentiable manifold 911 and a stationary, 
axisymmetric Lorentzian metric g with signature (1, 3), which in Boyer-Lindquist coordinates (t, r, 9, ip) [2], where 
t e IR, r € IR>o, 9 € [0, tt], and p € [0, 27r), is given by 

g = ^ (dt — asin^ (9) dtp) 0 {dt — asin^ (9) dtp) 

^ . 2 ( 1 ) 
— dtp — o dt^ 0 ^ [r^ -I- a^] dtp — a dt^ — ^ dr tg) dr — S d0 tg) d0 . 


The horizon function is defined by A = A(r) := (r — r+)(r — r_) = r^ — 2Mr + a^, r± := M ± V— a? 
denote the event and the Cauchy horizon, respectively, M is the mass and aM the angular momentum of the 
black hole with 0 < a < M, and S = S(r, 0) := r^ -I- cos^ (9). We describe Kerr geometry in terms of a 
local Newman-Penrose null tetrad frame that is adapted to the principal null geodesics, i.e., the tetrad coincides 
with the two principal null directions of the Weyl tensor. In this so-called Kinnersley frame [24], since Kerr 
geometry is algebraically special and of Petrov type D, we are presented with the computational advantage that 
the four spin coefficients K,a,X, and v vanish and only one Weyl scalar, namely 'I' 2 , is non-zero. Accordingly, 
the congruences formed by the two principal null directions must be geodesic and shear-free [26]. We construct 
the Kinnersley tetrad directly from the tangent vectors of the principal null geodesics [8] 


dt -I- 

“i— = -j-A 

dx A 


^ = ±E, 

dx 



d(p 

dx 


- A^’ 


( 2 ) 


where x is an affine parameter and E denotes a constant. Hence, aligning the real-valued Newman-Penrose 
vectors I and n with the principal null directions and choosing conjugate-complex Newman-Penrose vectors m 
and m in such a way that they satisfy the Newman-Penrose conditions (A1)-(A3) yields 




^ [r^ -\- a^]dt -\- ^dr a 


n = 


^-^^(^[r^ + a^]d,-Adr + ad^) 


m = 


'J2 {r -\- ia cos (0)) 


(ia sin {9) dt + dg +i esc {9) d^f 


( 3 ) 


m = — 


(r — iacos (9)) 


(ia sin {9) dt — dg + i esc (9) 


with the signum function 


sign(A) := 



for A > 0 
for A < 0 . 


Note that the use of the absolute value and the signum function allows for a unihed representation of the frame 
both outside and inside the event and the Cauchy horizon, respectively. For the calculation of the corresponding 
spin coefficients, i.e., for solving the first Maurer-Cartan equation of structure (A5), we require the dual co-tetrad 
of (3) given by 

I = sign(A) (dt — dr — a sin^ (9) dp 
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m = 


'J2 (r + ia cos (9)) 


m = — 



1 ( 

\/2 (r 

— iacos (0)) ^ 


^ia sin (0) dt — S d0 — i [r^ + a^] sin ( 

^ia sin (9) di + E d0 — i [r^ + a^] sin (9) d(p^ . 


We introduce further computational simplifications by applying a class III local Lorentz transformation (A7) 
with parameters of the form 


c = 


^ and = 
2E 


7e 

r — ia cos (9) 


to the Kinnersley tetrad (3). This transformation leads to the so-called Carter (symmetric) frame [4], which acts 
under the composition of the discrete time and azimuthal angle reversal isometries 11—>■ — t and (p i—>■ —(p as 

V I—>■ —sign(A) n', n' i—> —sign(A) V , m! i—> fn' , m' i—> ml 

and, thus, has only six independent spin coefficients 


k' = —v' , tt' = —t' , a' = —jS' , a' = sign(A) A', p! = sign(A) q' , e' = sign(A) 7'. 


We obtain 




(4) 


m' = —i= (ia sin (0) dt + dg +i esc (0) d^) 
v 2 E 


(ia sin (0) dt — de + \ esc (0) 




The dual co-tetrad reads 


I' = t/!^ sign(A) (dt— ^ dr — a sin^ (0) d(/? ) 
\ Z 2-i \ ZA y 


n' = ( dt -I- ^ dr — asin^ (0) d^?) 

V z 2-i \ za y 




^ia sin (0) dt — E d0 — i [r^ -I- a^] sin (0) d^?^ 


m' = — (i® sin (0) dt -I- E d0 — i [r^ -f a^] sin (0) dp 


(5) 


We point out that Boyer-Lindquist coordinates have singularities at the event and the Cauchy horizon. The light 
cone of an observer that, for example, approaches the event horizon from outside the black hole closes up and 
becomes degenerate (see Figure 1). Moreover, space and time reverse their roles in between the event and the 
Cauchy horizon. All this precludes a study of the propagation of Dirac waves across the horizons. In order to 
have a consistent description of their overall dynamics, we use advanced Eddington-Finkelstein-type coordinates 
(see [13, 14] for the original Eddington-Finkelstein null coordinates). This analytic extension covers the black 
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hole exterior and interior and allows for regular transitions of Dirac waves across the horizons. Furthermore, it 
possesses a proper coordinate time, which is required for the Hamiltonian formulation of the Dirac equation and 
the corresponding Cauchy problem. The advanced Eddington-Finkelstein-type coordinates are constructed as 
follows. By means of the tangent vectors (2), we derive two relations between the time and the radial coordinate 
and two relations between the azimuthal angle and the radial coordinate along the principal null geodesics 


dt _ 

dr A 


t = zL 


A 


■ dr + c± = ±r* + c± 


( 6 ) 


where 


d(/? ^ a 

dr A 


(p = ± 


® 1 / , o 1 

— dr + c+ = ±- In 

A r+ — r_ 


r — r+ 
r — r_ 


■c± 


r+ — r_ 


In |r — r+| — 


In |r — r_ I 


r+ — r_ 


is the Regge-Wheeler coordinate and c ±, c± are constants of integration, motivating a transformation to coordi¬ 
nates that are adapted to ingoing null geodesics 


M X IR>o X [0, tt] X [0, 27r) — >■ M x IR>o x [0, tt] x [0, 27r), 


(t,r,6»,(/?) (T,r,6>,^) 


with 


T := t + Ti, — r = t- 


r+ — r_ 


In |r — r_|_| — 


r+ — r_ 


In |r — r_ I 


(j):=ip + 



In 


r — r+ 
r — r- 


(7) 


In these so-called advanced Eddington-Finkelstein-type coordinates, light cones are non-degenerate at the hori¬ 
zons (see Figure 2). Approaching again the event horizon from outside the black hole, the light cone tips over 
until its outgoing future light cone is aligned with the horizon, which indicates its trapping characteristic. Be¬ 
sides, ingoing light rays are represented simply by straight lines. This can be directly seen from the two relations 
between the advanced Eddington-Finkelstein-type time and radial coordinates for ingoing and outgoing null 
geodesics 


— = —1 and 

dr I in 


dr 4M r 

IW = 1 ”1 -) 

dr lout A 


which are obtained by inserting (7) into (6). The metric (1) represented in these coordinates becomes 


r — asin^ (9) d(/)]^ 

/ 2Mr\ 

— M -I-j (dr — a sin^ (9) d<j>) 0 (dr — a sin^ (9) dcj)) —T,d9^d9 — Y, sin^ (9) d(j>^d(j). 

Considering the restriction to constant-r hypersurfaces reveals that these are space-like and that t is a proper 
coordinate time. In the Carter-Penrose diagrams shown in Figure 3, we depict the constant-t and constant- 
r hypersurfaces in Boyer-Lindquist coordinates and the constant-r and constant-r hypersurfaces in advanced 
Eddington-Finkelstein-type coordinates for Kerr geometry. While the constant-t hypersurfaces become time-like 
inside the black hole in region II, the constant-r hypersurfaces are always space-like and smoothly continued 
across the horizons. The Carter tetrad (4) and its dual (5) expressed in advanced Eddington-Finkelstein-type 
coordinates read 

l' = ([^ + 4Mr] dr + ^dr + 2a 

V 2 E |A| 


S= 1- 


/ 2Mr\ 2Mr / 

( 1-^ j dr 0 dr-^ (^ [dr — a sin^ (0) d^] 0 dr -I- dr 0 [d 


n = 




6 



FIG. 1: Causal structure of the non-extreme Kerr geometry in Boyer-Lindquist coordinates. A projection onto the (t, r)- 
plane, where every point is a 2-sphere, is presented. The real-valned Newman-Penrose null vectors I and n, which point 
along the principal nnll directions, form the light cones. The light cone of an observer approaching the event horizon from 
outside the black hole (r \ r+) closes up and becomes degenerate. In contrast, it opens up when the observer approaches 
the event horizon from inside the black hole (r Z' r+). This stems from the fact that the roles of space and time are 
reversed in the black hole interior region II. When r —>■ oo, the light cone becomes a 45°-Minkowski light cone because 
the spacetime is asymptotically flat. Note that all figures are restricted to regions I and II in order to avoid the issues 
that arise when one considers the ring singularity at (r = 0, 0 = 7r/2) and the maximum analytic extension. 



FIG. 2: Gausal structure of the non-extreme Kerr geometry in advanced Eddington-Finkelstein-type coordinates. A 
projection onto the (r, r)-plane is presented. Ingoing light rays are straight lines pointing in the n-direction. The light 
cone of an observer moving toward the event horizon from outside the black hole tips over until - after having crossed 
the event horizon - its future light cone is completely in the black hole interior. This shows the trapping characteristic of 
event horizons. 
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FIG. 3: Carter-Penrose diagrams for the non-extreme Kerr geometry in Boyer-Lindquist coordinates (left) and advanced 
Eddington-Finkelstein-type coordinates (right). The blue lines represent constant-t hypersurfaces, the red lines constant-r 
hypersurfaces, and the green lines constant-r hypersurfaces. The constant-t and constant-r hypersurfaces are restricted 
to either the exterior or the interior of the black hole. Their nature changes across the event horizon, i.e., space-like 
hypersurfaces become time-like and vice versa. The constant-r hypersurfaces (cut-off at the Cauchy horizon) are space¬ 
like outside and inside the black hole and smooth across the event horizon. 


m = 




(ia sin (0) dr + de +i esc (9) d^) 


and 


m' = — ^ — (ia sin {9) dr — de + i esc {9) d^) 
V 2 S 


i' = y^sign(A)(dr + 




dr — asin^ (9) dej) 




n' = \ (dr -I- dr — a sin^ (0) d^) 




^ia sin [9) [dr -|- dr] — S d0 — i [r^ -|- o^] sin (0) di))^ 


l^ia sin {9) [dr -|- dr] -I- S d0 — i [r^ -|- a^] sin {9) d())^ . 


Since the real-valued Newman-Penrose vector I' is still singular at the horizons, we apply a class III local Lorentz 
transformation (A7) with parameters 

^ ^ and Ip = Q . 

r+ 































This leads to the regular Carter tetrad 


I” = — ([A + 4Mr] + A + 2a a^) 

\/2 S r+ 




: (la sin (0) aT- + a^ + i esc (0) a^) 


m" = - - - (ia sin (0) a,- — a^ + i esc (0) a^) 

v2 S 


and to the dual co-tetrad 


l” = ^ — dr -I- 1 — ^ dr — a sin^ (0) d(/) 


n" = (dr -I- dr — a sin^ (0) d(/)) 


m = —- 


: ^ia sin (0) [dr -I- dr] — S d0 — i [r^ -|- a^] sin (0) d(/)^ 

= ^ia sin (0) [dr -|- dr] -|- S d0 — i [r^ -|- a^] sin (0) d(/)^ . 


Substituting the latter into the first Maurer-Cartan equation of structure (A5), we obtain regular spin coefficients 
for the non-extreme Kerr geometry in advanced Eddington-Finkelstein-type coordinates 


-P” = - (^) - sin(0)) , 


ia sin (0) „ r+ „ A 

-\/2 E (r — ia cos (0)) ’ ^ -\/2 E (r — ia cos (0)) ’ ^ -\/2E r+(r — iacos (0)) 


2^/2 -^E (r — ia cos (0)) ’ 


r^ — a^ — 2ia cos (0) (r — M) 
2^/2 -^/e r+ (r — ia cos (0)) 


In the next section, we show the separability of the massive Dirac equation in this horizon-penetrating frame, 
perform asymptotic analyses of the solution of the resulting radial system at space-like infinity, at the event 
horizon, and at the Cauchy horizon, and discuss the eigenfunctions as well as the eigenvalues of the resulting 
angular system. 

III. THE MASSIVE DIRAC EQUATION IN THE ANALYTICALLY EXTENDED KERR 

GEOMETRY 

A. Mode Ansatz and Separability 


Substituting the regular Carter tetrad (8) as well as the associated spin coefficients (9) into the massive Dirac 
equation (B5) and employing Chandrasekhar’s mode ansatz (see, e.g., [8]) 

^i(r, r, 0, (/>) = (r — iacos (0)) ^(r, 0) 

'^i{T,r,e,<j)) = (r -|- iacos (0)) ^i{r,0 ), 


( 10 ) 



9 


where w € IR is the frequency, k € Z +1/2 the wave number, and z G {1,2}, we obtain the first-order PDE system 

(^ Adr + r — M — iLO {A + AMr) — 2iak)j1fi + (dg + 2~^ cot (O) — auj sin (9) — k esc (d))^ 

= V^i/x* (r — ia cos (9)) 

r+(9r + iw) — {dg + 2~^ cot (0) -I- aw sin (9) + k esc {9))jifi = —•\/2i^*(r — iacos (9)) ^2 

{Adr + r — M — \u! {A + AMr) — 2iak) ^2 — {dg + 2~^ cot (9) + aw sin {9) + k esc (d)) 

= •\/2i^*(r -I- iacos (d))^ 


r+(9j. -|- iw) + {^dg + 2 ^ cot (9) — aw sin (9) — k esc ( 6 *)) ^2 = —V2i^^(^r -b iacos {9))jfi . 

This system is separable by means of the product ansatz 

Mi{r,9) = <^+{r)^+{9) 

J>f2{r,9) = ^-{r)^_{9) 

^i(r,0)=^_(r)b7+(d) 

^2ir,9)=M+{r)^.{9), 

leading to a radial ODE system 

{Adr + r — M — iu! {A + AMr) — 2iak)^g- = r+(^i /3 + •\/2i/i*r)^_ 
r+{^dr + iw)^_ = (^2/4 — •\/2i^*r)^+ 
and to an angular ODE system 

{dg + 2~^ cot ( 0 ) — aw sin ( 0 ) — fcese ( 0 ))l?h = —(^ 1/4 — •\/ 2 yx*acos ( 0 ))=^ 

{dg + 2~^ cot ( 0 ) + aw sin ( 0 ) + k esc ( 0 )) = (C2/3 + V2^i,a cos ( 0 )) .S^- 


( 11 ) 


( 12 ) 


with the constants of separation j G {1,2, 3,4}. Erom the radial system, it can be directly seen that the 
identifications = ^3 and ^2 = ^4 have to hold, whereas from the angular system, we obtain the identiheations 
Cl = C4 and C2 = Cs- Thus, defining ^ := ^3 = ^2 = Ca = ^4) (H) and ( 12 ) reduce to 


and 


{Adr + r — M — iuj {A + AMr) — 2iak)^+ = r+ (C -b •\/2i^*r).^_ 
rjr{dr + iw)^_ = (C — •\/ 2 i/x*z’).^+ 

{dg -b 2~^ cot (0) — aw sin (0) — fcese (0)) = — (C — •\/ 2 ^*acos {9))3^+ 
{dg + 2~^ cot ( 0 ) -b aw sin ( 0 ) -b k esc ( 0 )) = (C + ■\/ 2 /x*a cos ( 0 )) , 


(13) 


(14) 


respectively. The radial system (13) can be transformed into a more symmetric form by means of the functions 
= \/|A| and = r+ , yielding 


(a dr — iw (A -b AMr) — 2iafc).^+ = \/|A| (C -b \/2i^i,r)3S- 
A{dr -b iw).^_ = sign(A) ^|A| (C — . 


(15) 


For the following study of the asymptotics and the decay of the radial solutions at infinity, the event horizon, 
and the Cauchy horizon, it is convenient to rewrite (15) in the matrix representation 


dr^ = U{r)^ , 

where ^ := {^+,^-)^ and 

, . _ j_ / i (w(A -b 4Mr) -b 2ak) a/|A| (C -b v^iAi*?") \ 
^ ■ A V sign(A) a/|A| ($ - V^i9-*r) -iwA J 


(16) 
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This representation also allows for a straightforward determination of the singular points of the radial system [9]. 
We find that U has singularities of rank /i = 0 at r = r±. Therefore, even though both the coordinate system 
and the tetrad frame are regular at the event and the Cauchy horizon, these are regular singular points of (15). 
Thus, the regions outside the event horizon r_|_ < r, between the event and the Cauchy horizon r_ < r < r+, 
and inside the Cauchy horizon r < r- have to be considered separately. Note that this is a particularity of the 
mode ansatz (10). Finally, the matrix representation of the angular system (14), with := , is given 

by 


cos {9) 

dg + 2“^ cot (0) + aw sin {9) + kcsc{9) 


-dg - 2- 


cot {9) + aoj sin {9) + kcsc{9) \ _ c 

cos (9) J 


(17) 


B. Asymptotic Analysis of the Radial Solution at Infinity 

Following the approach of [17], we derive the asymptotic solution of the radial system (16) for r —)■ oo and examine 
its decay. This solution is required (in addition to the asymptotics at the event and the Cauchy horizon, which 
are determined in the subsequent subsections), e.g., for a description of the scattering process of Dirac waves 
by the gravitational field of a Kerr black hole (cf. Section IV) and for the construction of a functional analytic 
integral representation of the Dirac propagator [18]. We begin by expressing (16) in terms of the Regge-Wheeler 
coordinate 


= T(r^)M, (18) 

where T{ri,) := A/(r^ + a^) U{r-i,). Then, we represent this radial system by the diagonal matrix S := D~^TD = 
diag(Ai, A 2 ), with A 1/2 being the eigenvalues of T, and by the corresponding diagonalization matrix D 

dr,{D-^^) = [S-D-^ {dr,D)] . 


Using the ansatz 


.^(r*) = D(r*) 


ei‘^+('-*)fi(r,,) 

g-i0_(r*)^2(r*) 


we obtain an ODE system for f := (fi, f 2 )^ 

drj= [S-W-^D-^ {dr,D)W -W-^dr,W]f, 

where W := diag(e''^+,e““^“). The functions (j)± are fixed by requiring that S = W~^dr^W, i.e., 

dr^4>+ = —iAi and dr^4>- = iA 2 , (19) 


yielding 

drJ=-W-^D-^ {dr,D)W^. 

Lemma III.l. Every nontrivial solution ^ of (18) is asymptotically os r —)■ 00 of the form 


~ ~ / e^4>+(r*) fW \ 

^{r^) = + E^{r^) = ^( 2 ) ) + L^oo(?'*) 


with the asymptotic diagonalization matrix 

cosh (O) sinh (O) 
sinh (D) cosh (fl) 




1 j cosh (O) + i sinh (fl) sinh (fl) + i cosh (D) 
[ y/2 I sinh (D) + i cosh (fl) cosh (O) + i sinh (D) 


for > 2p,l 


for oj^ < 2pl , 


( 20 ) 


( 21 ) 


(22) 
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where 


Q. := I 


In 


- In 


OJ - v^/i* \ 

for > 2p,‘: 

oj y 

-uj\ 

for < 2p 

■\/2/i* + bJ J 


the asymptotic phases 


(/)±(r*) ~ sign(w) x <( 


— — 2^2 _|_ 2M ± 0 ;- * In (r*) for uP > 2/rJ 

V - 2^5 / 


— a;2 r* + 2M ±uj -) In (r*) for uj‘^<2plI, 




and an error E^o with polynomial decay 


||£^oo(?'*)|| = \\p?{r*) - ^oo{rP\\ < — 

for a suitable constant a G M>o • The asymptotics of the function f for large values of r is given by 

^T 


foo := = const. 7^0 


with an error Ef that has polynomial decay 

\\E^{rP\\ = \WP-U\<^ 

for a suitable constant b G K>o- 

Proof. The matrix T defined below (18) converges for r —>■ oo to 


Too '■= lim T = i 

r* —>-oo 




(23) 


(24) 


(25) 


(26) 


Further, it has a regular expansion in powers of 1/r*, i.e., T = Too+0(l/r*). As a consequence, both the diagonal 
matrix S and the diagonalization matrix D also have regular expansions in powers of 1/r*. The eigenvalues of 
(26) read 


Ai /2 ^ sign(a;) x 


=F \\/wP — 2/x* G C 
=F-\/2/x* — up G M 


for bP > 2p,l 
for < 2/iJ . 


The associated matrix Doo '■= hrur^-yoo D is given by (22) with arguments (23). This can be easily verified by 
direct calculation. With the asymptotic first-order eigenvalues of T 


■ ^ 2iM / 

T _ 2^2 ^ ^ 

r* \ 




a/w 2 -2/r2 


Ai /2 - sign(w) X 


T^2p.l - + — (iwi 


^2iP-uj^ 


for > 2/rJ 


for < 2/iJ , 


we can directly solve (19) by simple integration and obtain (24) for the asymptotic phases. Next, as the Hilbert- 
Schmidt norms of D~^ and dr^D are bounded from above by 


\\D ^||hs<c and 


||9r*T>|lHS < 


d 


r 


2 
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for r* sufficiently close to infinity, where both c and d denote positive constants, we estimate the C^-norm of 
(20) by 


||5.J||<2p-l||HS-||9.*l?|lHS-|lfl!< 


(27) 


with ||bb||HS = IIbb ^||hs = V^- Employing the triangle and the Cauchy-Schwarz inequality, we derive the 
following inequality 

I ia..(f,f)| l(f,a.j) + (9,j,f)i |(f,9.j)l + l(a,j,f)| |(f,5,j)i ||fll-lia,j|| 


l^r-JIfl" 2 ||j^|j 2 I 

Consequently, using (27), we get 


< 




< 


= l|a.J||. 


(28) 


Note that ||f|| ^ 0 because ^ has to be nontrivial. Integrating (28) over the Regge-Wheeler coordinate from tq 
to r* and applying the triangle inequality for integrals gives for all 0 < ro < r* 

r*d< 

„/2 

\J ro I ro Jro ’ * 

and, hence. 


dr' In IlfII dr; 


< / \dr' In llfll dr* < 2cd 




2cd 

In llfll 


< - 

— 



K 


Then, since 0 < 2cd/r'^'\° < oo, there exists a constant fV > 0 such that 

^ < llfll <A^- 

Substituting this into (27), we find for sufficiently large r* 

\\drJ\\<^,, 


(29) 


(30) 


where b := 2cdN. This implies that f is integrable and, according to (29), has a finite, non-zero limit foo := 
limr*_>oo f(''’*) 7^ 0- Finally, integrating (30) from r* to oo and again using the triangle inequality for integrals, 
we obtain the error estimate 


Pfll = llf-foo|| = 


pOO pOO 7 

/ 5^'fdr; < / lia^'flldr; < — . 

t/ r* 1/ r* 


(31) 


The l/r*-decay of the error Eao in (25) follows directly from (31). 


C. Asymptotic Analysis of the Radial Solution at the Event Horizon 

Using the solution ansatz 

V fl 2 (r.) J ’ 

( ~l“) '2 

with the angular velocity of the event horizon := a/ (2Mr+), in (18) yields an ODE system for g := (fli, 02 ) 


dr,9 = - 


j,2 _|_ q 2 


2k r - a 


1 0 
0 0 


Va 


/ Va{u; + 2kn^+l) (^2^*^ - i^) \ 


Ve2iM<->* (y2^*r + iC) 


^/A 


UJ 




(32) 


0- 
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In the limit r r+j the right-hand side of this system vanishes and, thus, we obtain the asymptotic solution 
0r+ := 0 = COnSt. 

Lemma III.2. Every nontrivial solution ^ of (18) is asymptotically as r \r+ of the form 


^(r*) = ^r+ (r*) -I- Er^ (r*) 


^2iL+kn^+J)r, ( 1 ) 

- ^ ' }^r+ 

9? 


with 

0r+ := (01+^0^+^^ = const. 7 ^ 0 
and an error Er_^ with exponential decay 




+ £^r+(c*) 


(33) 


for r sufficiently close to r+ and suitable constants p,q G K>o- 

Proof. For r \ r+, as r ~ r+ -|- with q := (r+ — r_)/(4Mr+) € ffi>o, the right-hand side of (32) is of the 
order ©(e'^’’*). Hence, there exists a constant p' G M>o such that for r* sufSciently close to —oo 


||9.,0|| <pV'-*||fl|i. 


(34) 


Similar to the proof of the previous subsection, we find for all r* < tq 


ln|l 0 | 


< ?- e^’’* I 
- q h* 


and since e^'’* is positive, there is a constant N' > 0 such that the C^-norm of g is bounded 

1 


N' 


<ll0ll<iV'. 


Substituting (35) into (34) yields 


\\drM\<Pe'^^% 


(35) 

(36) 


where p := p'N'. All this implies that g is integrable and has a finite, non-zero limit for r* —>■ 
(36) from —oo to tq and applying the triangle inequality for integrals, we obtain 


\\Es\\ = Il0-0r+|| 



< 


Pwgll dr; < pe'?''* 


—oo. Integrating 


(37) 


which proves the exponential decay of the error E^. The exponential decay of the error Er_f_ follows directly from 
(37). ■ 


D. Asymptotic Analysis of the Radial Solution at the Cauchy Horizon 

We begin with the solution ansatz 

V Ur.) ) ’ 

where := a/(2Mr_) is the angular velocity of the Cauchy horizon, and employ it in (18). This leads to an 

ODE system for 1) := (()i, [} 2 )^ 


= 


r2 _|_ fj2 


-2k{2Mil^-lr-a)(l o)+^ 


^ ^ (a; + 


e-2 i(v^p*r-iC) ^ 

vlA| w / 


If. 
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In the limit r r_, the square bracket on the right-hand side vanishes, resulting in the asymptotic solution 
^ = COnst. 

Lemma III.3. Every nontrivial solution ^ of (18) is asymptotically as r \r-. of the form 


^(r*) = (r*) -I- Er_ (r*) = 


2 i 


(-) 


[,( 1 ) \ 


ll?_> 


with 


t)r_ := = COnst. 7^ 0 

and an error Er_ with exponential decay 

\\Er_{r.)\\<ue-^^* 


for r sufficiently close to r_ and suitable constants u,v G M>o. 

Proof. The proof of this lemma is similar to the proof of Lemma III.2. ■ 


E. Eigenfunctions and Eigenvalues of the Angular System 

The angular system (17) in its decoupled second-order form is known as the massive Chandrasekhar-Page equation 
[8]. In the limit a \ 0, the solutions of this equation reduce to the spin-weighted spherical harmonics for the 
spin-1/2 case [21], For non-zero angular momenta, the solutions are usually referred to as the spin-1/2 spheroidal 
harmonics. For a good introduction and a compilation of some properties of these functions, we refer the reader 
to the recent paper [10]. In the present work, it is only of importance that the matrix-valued differential operator 
on the left-hand side of (17) has a spectral decomposition with discrete, non-degenerate eigenvalues and smooth 
eigenfunctions, which was proven in [15, 17]. In the following, we state the result. 

Proposition III.4. Eor any w € IR and k G + 1/2, the differential operator in (17) has a complete set of 
orthonormal eigenfunctions m L^((0, tt), sin (0) d6*) . The corresponding eigenvalues fn real-valued 

and non-degenerate and can, thus, be ordered as fn < Cn-i-i- Moreover, the eigenfunctions are bounded and 
smooth away from the poles. Both the eigenfunctions and the eigenvalues depend smoothly on uj. 


IV. SCATTERING OF MASSIVE DIRAC WAVES BY THE GRAVITATIONAL FIELD OF A KERR 

BLACK HOLE 


In this section, we study the scattering of massive Dirac waves by the gravitational field of a Kerr black hole 
from the point of view of an observer described by horizon-penetrating advanced Eddington-Finkelstein-type 
coordinates. To this end, we consider Dirac waves that emerge from space-like infinity and propagate toward 
the black hole’s event horizon. We compute the net current at infinity as well as at the horizon, and derive 
a conservation law for the reflexion and transmission coefficients. Starting with the asymptotic radial solution 
at infinity (21) and at the event horizon (33), we impose boundary conditions specifying an incident wave of 
unit amplitude at infinity, a reflected wave of amplitude at infinity, and a transmitted wave of amplitude 

at the event horizon. The asymptotic ingoing and outgoing radial solutions adapted to these boundary 
conditions read 


and 


~ f A p^'P+P*) \ 

^Scat.(r^oo)cs j 

^Scat.(?’ \ ?'+) ~ ^ ^ . 


(38) 


(39) 


We point out that the boundary conditions are chosen in conformity with the physical requirement that no 
waves can emerge from the event horizon. Besides, only the asymptotic solution branch of (21) with > 2pf 
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is considered because free particles at infinity must have energies that exceed - or at least equal - their rest 
energies. Next, assuming the normalization condition ||^(0)|p + ||.^(0)|p = 1 for the angular eigenfunctions 
(cf. Proposition III.4), the radial Dirac current (for dehnitions and notations see Appendix B) yields 


r = ^ (sign(A) ||^+f - 

with the radial Infeld-van der Waerden symbol 


AB 


r m^\ _ 1 

fn^ 


A 0 \ 

0 -^IJab 


As this current has jump discontinuities at the event and the Cauchy horizon (see the final paragraph of Section 
III A), the regions separated by the horizons have to be evaluated independently. For the scattering problem at 
hand, however, it suffices to consider the exterior region of the black hole r+ < r, where the radial Dirac current 
reads 


r 




(40) 


From the radial system (15) and the corresponding conjugated complex system, we obtain the relation 


^ ^ = const. 


by simple algebraic manipulations. Substituting this into (40), we derive the conserved net current of particles 


dtN = 




J'' ^|det(g)| d9 dcj) 




^ = const., 


where N is the total number of particles and ■\/|det(g)| = Esin(0). Defining the reflexion and transmission 
coefficients 








and using (38) and (39), the net current at infinity and at the event horizon becomes 


dtN 


1^- 


r+ 


(41) 


and 

Att 

dtN\r\r+ = - ) (42) 

r_|_ 

respectively. It follows from the latter equation that the net current across the event horizon is always positive. 
Furthermore, with the constancy of the net current and (41) and (42), we deduce 

= 1 ) 

which proves that superradiance cannot occur because the reflexion coefficient is always less than unity. These 
results are in agreement with those found in the accordant analysis of this scattering problem employing Boyer- 
Lindquist coordinates (see [8] and references therein). 
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Appendix A: The Newman-Penrose Formalism 


Let (911, g) be a Lorentzian 4-inanifold endowed with the unique, torsion-free Levi-Civita connection a; and dual 
basis (e^) and (e^), g G {0,1,2,3}, on sections of the tangent and cotangent bundles TDJl and T*91t, respectively. 
Further, let FDJl and F*91l be local (orthonormal or null) frame bundles on 911. The dual basis on sections of FdJl 
and F*dJl each consist of four vectors denoted by (e(a)) and (e(“)), a G {0,1, 2, 3}. In terms of the original basis 
vectors, these can be written as and 6 ^“^ = where is a linear map from r91l to F91t, 

namely a (4 x 4)-matrix, and is its inverse. In the Newman-Penrose formalism [25], the local basis vectors 
are given by two real-valued null vectors, I = and n = as well as by a conjugate-complex 

pair of null vectors, m = 6 ^ 2 ) = —and fn = These have to satisfy the null conditions 

ll = nn = mm = fnrn = 0, (Al) 

the orthogonality conditions 


I ■ ra = I ra = n ■ m = n rfi = 0, 

and the cross-normalization conditions (which depend on the signature convention) 

I ■ n = —m • m = I. 


(A2) 

(A3) 


The local metric in the Newman-Penrose formalism is non-degenerate and constant. It reads 
V = ® = I ® n + n®l — m®rn — ra® m . 

Since the Levi-Civita connection is torsion-free, the first Maurer-Cartan equation of structure simplihes to 

de(“) 


(b)(c)' 


(A4) 


where the Ricci rotation coefficients are related to the connection 


via 




' p(c) _ i®r\p^ -4-p 

(b)(c)® “ (6) + fi/i 6 (6)^ p- 


In the Newman-Penrose formalism, (A4) becomes 

dl = 2?ft(e) nA I — 2nA ?ft(Kfn) — 21 A 3?([t — a — P] fnj + 2i9(p) m Arn 


dn = 23 ^( 7 ) n Al — 2n A 3?([a -I- /3 — tt] m) -|- 2i A 3?(i^m) -|- 2i9(^) m Am 


(A5) 


dm = dm = {tt + t) n A I + (2i SI(e) — g)nAm — anAm+(^g + 2i 9 ( 7 )) I A m 
+ XI Afn — {a — P)m Am 


with the so-called spin coefficients 

= 9'(2)(0)(0) 

9 = 7(2)(0)(3) 

^ - 2(7(1)(0)(0) +7(2)(3)(0)) 


= 7(2)(0)(2) 

9 = 7(i)(3)(2) 

+ 

o' 

.-HiC'I 

II 


^ = 9'(1)(3)(3) 

^ = 7(2)(0)(1) 

« - K7(1)(0)(3) +7(2)(3)(3)) 

(A6) 

= 7(1)(3)(1) 

^ = 7(1)(3)(0) 

P = M7(1)(0)(2) +7(2)(3)(2)) 



and the real and imaginary parts 3?(-) and SI(-), respectively. The transformations employed in this study are 
elements of the 2-parameter subgroup of local Lorentz transformations known as class III or spin-boost Lorentz 
transformations. These renormalize the real-valued Newman-Penrose vectors I and n, but leave their directions 
unchanged, and rotate the conjugate-complex pair m and fn by an angle ip in the (m,Tn)-plane [8] 

I i-A l' = , n PA n' = n , m pa m! = F'^m , m 1 —> m' = fn , (A-7) 

where c G IR\{0} and G M are functions depending on the spacetime coordinates (x^). There are various 
aspects of the Newman-Penrose formalism that are not discussed in this appendix such as the different classes of 
local Lorentz transformations or the Weyl scalars and their algebraic classification. They can be found elsewhere 
in the literature. The interested reader may be referred to [26, 28]. 
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Appendix B: The General Relativistic Dirac Equation 

The general relativistic, massive Dirac equation without an external potential is given by [19, 32] 

+ im)'I'(a:'^) = 0 , 

where the 7 ^ are the general relativistic Dirac matrices, which satisfy the anticommutator relations { 7 ^, 7 “^} = 
2 ( 7 ^'^idc 4 , is a Dirac 4-spinor on sections Sx^ — C"* of the spin bundle SVJl = 911 x C"* on 911, V denotes 

the metric connection on S'911, and m is the fermion rest mass. Using the chiral 2-spinor representation of the 
Dirac 4-spinors and matrices [20] 


4- = 


{Q 


0 


AB 


-fiAB 


^€{1,2}, Be{i,2}., 


with the 2-component spinors and Qg and the Hermitian (2 x 2)-lnfeld-van der Waerden symbols we 


iB 

obtain the following 2-spinor form of the Dirac equation 


^ AB^ +iM*QB — 0 


abQ T b — ^ ’ 


(Bl) 


where := m(\/2 and l^ote that dotted indices are subjected to conjugated complex 

transformations. Next, let and k € {1,2}, be (dual) local Newman-Penrose bases for the Dirac 2- 
spinors. The associated local spinor components can be expressed by and = C(k)^^A with 

the original 2 -spinor components as well as the (2 x 2 )-matrix and its inverse C(fc)'^- In this 

representation, the metric connection reads 






(n)(fe)(/) 


f^{n) 


where and 


^(m) 


(n){k)(i) 


= V 


(m)(6) 


{n){d)(k){l) 


— f(°)(P) rr’', , rr'’ • p p 

— vxe e a a cr e„ e 


. p 


(b) (c) 

7(a)(6)(c) 


(B2) 


(B3) 


We point out that the 2-dimensional Levi-Civita symbol e acts as skew metric on C^. Furthermore, the Infeld-van 
der Waerden symbols yield 


(k)(i) 


P mi^ 
fn^ 


By means of (B2)-(B4), (A6), and the definitions := '■= P^‘^\^i := Q^‘^\ as well 

the general relativistic Dirac equation (Bl) in the Newman-Penrose formalism becomes 

+ e - g)^i + ijn^d^ Pit - 0)^2 = i^*^! 

{n^d^ + M - 7)^2 + -f /? - r)^i = 

{Pdfj, +e- q)% - {m^d^ Pif - a)^i = i^*^2 

[rPd^ + it - 7)^1 - (tn^ 9 ^ + P - t)% = • 


as 


:= -Q 


(B4) 

(i) 


(B5) 
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